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Abstract 

We study spintronics and spincaloritronics in topological insulators. We show spintron- 
ics effects in 2D topological insulator junctions and 3D topological insulators coupled 
to ferromagncts. We also investigate spin polarization on the surface of a topological 
insulator induced by a circularly polarized light and refraction at the junction between 
two topological insulators. As for spincaloritronics effects, we show transverse magnetic 
heat transport and thermoelectric transport in topological insulators. Finite-size effect 
in topological insulators and band structure engineering of interface states of topological 
insulators are also discussed. 
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1. Introduction 

Recent discovery of topological insulator (TI) offers a new state of matter topologically 
different from the conventional band insulator. Edge channels or 

surface states of the TI are topologically protected, immune to small perturbations which 
respect time reversal symmetry, and described by Dirac fermions at low energies. The 
surface Dirac fcrmion of TI has been predicted to show interesting phenomena, e.g 



sup erconducting proximity effect, in particular Majorana fcrmions[8 
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Since on the surface of the TI, electrons obey the (massless) Dirac equation, the spin 
and momentum of the electrons are coupled. Thus, one may expect novel spintronics 



and spincaloritronics effects on the surface of the TI. 
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One salient feature of the Dirac fermions is that the Zeeman field acts like a vector 
potential. Therefore, one can expect anomalous spin related property by magnetic fields 
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in TIs. We will see this feature in TIs attached to ferromagnets. The coupling between 
spin and momentum also gives novel coupling between spin transport and heat transport. 
This novel spin-heat coupling emerges not only as spincaloritronic effects, but also as 
thermoelectric effects. Namely, the coupling between charge transport and heat transport 
is also unconventional in TIs. 

In this paper, we study spintronics and spincaloritronics in topological insulators. 
We show spintronics effects in 2D topological insulator junctions and 3D topological 
insulators coupled to ferromagnets. We also investigate spin polarization on the surface 
of a topological insulator induced by a circularly polarized light and refraction at the 
junction between two topological insulators. As for spincaloritronics effects, we show 
transverse magnetic heat transport and thermoelectric transport in topological insulators. 
Finite-size effect in topological insulators and band structure engineering of interface 
states of topological insulators are also discussed. 



2. Spintronics 

In this section, we first consider spintronics effects in 2D TIs (or quantum spin Hall 
(QSH) systems ). Then, we study spintronics effects in 3D TIs. Transport properties of 



TIs have been nicely reviewed in Ref. 



24|. 



2.1. Spin rotation in normal metal/ quantum spin Hall junctions 

Let us commence with the scattering problem in normal metal/quantum spin Hall 
(N/QSH) junctions and show how the spin of an electron rotates upon scattering at the 
interface (see Fi g, [p . The effective 4-band model proposed for HgTe/CdTe quantum 
wells is given by [25| 

n = V (!) 

V h*(-k) j 

withh(k) = e(k)h x2 + d a (k)a a ,e(k) = C-D(k 2 x + k 2 y ),d a (k) = (Ak x , -Ak y , M (fc)) , and 
M(k) = M-B{kl+kD where we have used the basis order \H X +), \E X -), \Hi~)) 

("E" and "H" represent the electron and hole bands, respectively), and, A, B, C, D, and 
M are material parameters that depend on the quantum well geometry. e(k) is the 
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Figure 1: N/QSH junction and corresponding band structures (below). Dotted lines represent helical 
edge modes. 



averaged energy dispersion, while M(k) is the energy difference between the valence and 
conduction bands. A represents the strength of the band mixing due to the spin-orbit 
interaction. We note that % is equivalent to two copies of the massive Dirac Hamiltonian 
but with a fc-dependent mass M ffc).}^ In this system, the transition of electronic band 
structure occurs from a normal to an inverted type when the thickness of the quantum 
well is varied through a critical thickness. This corresponds to the sign change of the 
mass M of this system. 26|, |27j It should be noted that since h(k) ^ h*(—k) (which 



corresponds to the fact that chirality is different for different spins), spin orbit coupling 

has z component, which results in spin rotation effect. 

By solving scattering problem in this junction with the use of the Hamiltonian Eq.JT]), 

one can calculate the rotation angle of the spin of the electron within the plane a. [28 1 

Figure [2] shows the results for a in the plane of (9, M) for C = —0.08 eV in (a), C = —0.1 

eV in (b), and C = — 1 eV in (c) with the other parameters fixed as A = 4 eV-A, B = — 70 

cVA 2 and D = —50 eVA 2 . Here, 9 is the angle of incidence as shown in Fig. [T] These 

parameters are appropriate for the HgTe/CdTe quantum well. [25[ In Figure [2] (a), a 

sharp ridge in M < and 9 > region and its negative correspondence in M < and 

9 < are seen. This is in stark contrast to the case of usual insulators M > although 

a is still nonzero there. Note that the height of the ridge is as high as ~ n/2. With the 
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Figure 2: Spin rotation angle a as a function of injection angle and M, (a) C = —0.08 eV, (b) C = —0.1 
eV and (c) C = —1 eV. Note that a = n and a = —ir are equivalent and the continuous increase of a 
in the clockwise direction is separated into several sheets in (b) and (c). Note that a(—0) = —a(9) is 
satisfied. 

increase of \C\, we find a qualitatively different structure. Near the origin in Fig. [2](b), 
a reaches 7T, changes its sign, and winds by 4tt around the origin, while it does not in the 
region far away from the origin. The angle at which the dispersions of the helical edge 



modes cross the Fermi energy is given by Oc = ± sin 



AID 



We find that 



[Ak F VB 2 -D 2 _ 

a large magnitude of a in Fig. [5] appears around this angle, which means that the helical 
edge modes resonantly enhance the spin rotation. A similar problem in 3D TIs has been 
also studied in Ref. [29|. 

2.2. Magnetoresistance in a ferromagnet/ferromagnet junction on the surface of a topo- 
logical insulator 

Next, let us proceed to 3D TIs coupled to ferromagnets. We consider 2D ferromag- 
net/ferromagnet junctions which is abbreviated as F1/F2 in the following. We focus on 
charge transport at the Fermi level inside the bulk gap of the TI, which is described by 
the 2D massive Dirac Hamiltonian 

H _ / m z k x +m x - i(k y + m y ) \ 

\ k x + m x + i(k y + m y ) -m z I 

where m x ,m y and m z are exchange field and we set vf = h = 1. Here, x-axis points 
to the interface normal while y-axis is parallel to the interface. We choose the exchange 
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Figure 3: Tunneling conductance <r for m2 = ((a) and (c)), and mi = V0.9E ((b) and(d)). n-n 
junction in (a) and (b). p-n junction in (c) and (d). 



field in the Fl side as mi = (m x ,m y ,m z ) — mi (sin 9 cos tp, sin# sin ip, cos 9) while in the 
F2 side, we set m x = m y = and m z = m,2- 

By using the Landaucr formula, one can calculate the tunneling conductance of the 



junction. 



30( | Below, "n" and "p" mean that the Fermi level crosses the upper and the 



lower bands, respectively 

Figure [3] depicts the normalized tunneling conductance a in n-n junction for (a) 
m2 = and (b) m-i = \/0.9E where E is the Fermi energy. In Fig. [3] (a), the F2 is no 
longer ferromagnetic. However, the conductance strongly depends on the direction of the 
magnetization in the Fl. At 9 = or n, the mismatch of the wavefunctions between the 
two sides and that of the sizes of the Fermi surfaces suppress the conductance, because 
the energy E is near the bottom of the upper band in Fl while there is no gap in F2. 
At 9 = tt/2, on the other hand, the wavefunctions and the sizes of the Fermi surfaces 



are the same on both sides except the shift of Fermi surface in the momentum space due 
to the in-plane component of the magnetization. However, this misfit of the in-plane 
momentum between the two sides leads to a strong dependence of a on the in-plane 
rotation angle ip, which is not seen in the conventional magnetoresistance effect. Since 
k y is conserved due to the translational symmetry along y-axis, the positions of the Fermi 
surfaces strongly affect the charge transport: if exchange field points to x-axis, there is 
no evanescent wave. On the other hand, when exchange field is applied in y-direction, the 
Fermi surface moves to the k y direction and hence the overlap region between k y 's in the 
Fl and F2 is reduced. Therefore, the number of the evanescent modes increases and hence 
the conductance is strongly suppressed. Thus, we can obtain giant magnetoresistance in 
this system. In Fig. 121(b), the conductance is large at the parallel configuration (6 = 0) 
while it is small for antiparallel configuration (6 = it) . Tunneling conductance in p-n 
junction is shown for ni2 — in Fig|3] (c) and m,2 = y/0.9E in Fig. [3] (d). In Fig. [3] 
(d), in stark contrast to the conventional magnetoresistance effect, the conductance takes 
minimum at the parallel configuration (9 = 0) while it takes maximum near antiparallel 
configuration (6 = 7r). These 6 dependences can be understood by the overlap intcrgral 
of the wavefunctions on both sides. [3^ 

2.3. Magnetization dynamics on the surface of a topological insulator 

Up to here, we have considered TIs with static magnetizations. We next consider 



dynamical magnetizations. In Rcf. 



31 



|, the dynamics of magnetization coupled to the 
surface Dirac fermions of a 3D TI has been investigated by deriving the Landau-Lifshitz- 
Gilbert equation in the presence of charge current. It is shown that both the inverse 
spin-galvanic effect and the Gilbert damping coefficient are related to the two-dimensional 
diagonal conductivity of the Dirac fermion, while the Berry phase of the ferromagnetic 
moment is related to the Hall conductivity. The derived Landau-Lifshitz-Gilbcrt equa- 
tion has complicated torque terms, reflecting the spin structure of the surface Dirac 
fermion. Anomalous behaviors in various phenomena, e.g., the ferromagnetic resonance, 



31 



are predicted in terms of this Landau-Lifshitz-Gilbcrt equation in Ref. 

Since there is a one-to-one correspondence between the direction of momentum and 
that of the spin on the surface of the TI, an injected charge current can induce the 



magnetization, which is the so-called inverse spin-galvanic effect. 
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Figure 4: Time evolution of the magnetization vector. Ux,y,z is 3 components of the magnetization unit 
vector. Q is the frequency of the applied voltage. 

of the inverse spin-galvanic effect is a direct consequence of the fact that, on the surface of 
TI, the velocity operator is given by the Pauli matrices in spin space. The current-induced 
magnetization on the surface exerts torque on the magnetization of the ferromagnct when 
they are noncollinear. Therefore, by injecting an electric current on the surface of the TI, 
one obtains spin transfer torque between the surface Dirac fermion and the magnetization. 



In Ref. [321] , dynamics of the magnetization coupled to the surface Dirac fermions of 
a 3D TI has been studied by numerically solving the Landau-Lifshitz-Gilbert equation 
in the presence of charge current. The current-induced magnetization dynamics and the 
possibility of magnetization reversal have been discussed. Figure HI shows an example 
of current- induced magnetization reversal from n x = 1 to n x = — 1. Since the magne- 
tization dynamics is a direct manifestation of the inverse spin-galvanic effect, another 
ferromagnet is unnecessary to induce spin transfer torque in contrast to the conventional 
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setup composed of two ferromagnets.|3; 




Figure 5: The current on the surface of the TI as a function of time t induced by the rotating exchange 
field. (f> is the angle of the precession axis measured from the normal to the interface. 



In Rcf. 



34| . a current dynamics on the surface state of a 3D TI induced by an 



exchange field of an attached ferromagnet has been discussed. It is found that the 
magnetic precession generates a direct charge current when the precession axis is within 
the surface plane as shown in Fig. [5] This rectification effect is due to a quantum 
anomaly and is topologically protected. Intuitively, this phenomenon may be explained 
by an analogy to the quantum Hall effect: The exchange field perpendicular to the 
surface of TI plays a role of a magnetic field in the quantum Hall effect while that 
parallel to the surface corresponds to an applied inplane electric field. When the sign 
of the perpendicular exchange field changes, then that of the velocity of the parallel 
exchange field also changes. Therefore, simultaneous sign change of the magnetic field 
and the applied electric field does not change the sign of the Hall current: the rectification 
effect. In Ref. [35|, a similar rectification effect has been studied on an edge state of 



two-dimensional TIs induced by a circularly polarized light. 



2.4. The inverse Faraday effect on the surface of a topological insulator 

Optical effects of a TI have been studied theoretically in Refs 
the inverse Faraday effect 
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Due to 



43j, spin polarization can be induced on the 



surface of TI by a circularly polarized light. In Rcf. 



44j | . the spin polarization of the 



electrons on the surface of TIs under a circularly polarized light has been investigated. 



It is predicted that a light illumination can induce the out-of-plane component of the 
spin polarization as a result of the inverse Faraday effect. The magnitude of the spin 
polarization is proportional to the square of the lifetime r and (a z ) ~ 2 x 1CP 10 A~ 2 for 
typical parameters. 



2. 5. Refraction at the junction between two topological insulators 

Because the surface states of TIs are spin-filtered, they are expected to show novel 
charge transport which cannot be realized in other systems. We note that the linear 
dispersion of the Dirac cone on the surface of the TI is similar to photons, and the 
velocity of the Dirac cone on the surface of 3D TI depends on materials: about 4 x 10 5 m/s 



for Bi 2 Te 3 [52|, and 5 x 10 5 m/s for Bi 2 Se 3 



51 



Therefore, when two different TIs 



are attached, the refraction phenomenon similar to optics is expected at the junction. 
The effective Dirac Hamiltonian of the surface states on the xz-plane is represented as 
H = v(a x k)j,, where a arc the Pauli matrices, and v is the Fermi velocity. This 
Hamiltonian has the linear dispersion E = ±vk where k = |k| is the wavenumber. From 
this Hamiltonian we can easily see that for v > (v < 0), the spin on the Fermi surface 
on the upper cone has the clockwise (counterclockwise) direction in k space. Hence the 
sign of v corresponds to the spin chirality in k space. In Ref. j^ . a refraction problem 
is considered for a junction between the two TIs, which we call Til and TI2, with the 
incidence angle 8, and the transmission angle 9' (Fig. H3(a)). Let v\ and u 2 denote the 
velocities of the two TIs. 




(b) fefo=0.S 



(c) vJvi=1 .4 
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Figure 6: (a) Schematic of the refraction of the surface states at the junction between the two TIs, Til 
and TI2. (b)(c): Reflectance (red) and transmittancc (blue) for the ratios of the velocities of the two 
TIs: (b)v2/vi = 0.6 and (c)v 2 /vi = 1.4. 
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We assume v\ and V2 to be positive. After a straightforward calculation, the trans- 
mittance and the reflectance are given as 

sin 2 ^ cos cos 0' 

The results are plotted in Fig.[5£b)(c). Unlike optics, we have perfect transmission (T = 1 
and R = 0) for normal incidence (0 = 0) [5^ |. which reflects the prohibited backscattering 
on the surface of the TI, corresponding to the Klein paradox proposed in graphene 55^ . 

3. Spincaloritronics 

Here, we focus on transverse and longitudinal heat transports in TIs in Sees. 3.1 and 
3.2, respectively. 

3.1. Transverse magnetic heat transport on the surface of a topological insulator 

In Ref. [45|, Nernst-Ettingshausen and thermal Hall effects have been investigated 
on the surface of a TI on which a ferromagnet is attached. General expressions of the 
Peltier and the thermal Hall conductivities arc derived, which are reduced to simple forms 
at low temperatures. It is shown that the Peltier and the thermal Hall conductivities 
show non-monotonous dependence on temperature. At low temperature, they have linear 
dependence on temperature. From the behavior of the Peltier conductivity at low tem- 
perature, one can estimate the magnitude of the gap induced by time-reversal symmetry 
breaking. The Peltier conductivity can be used to map the Berry phase structure. 

3.2. Thermoelectric transport in topological insulators 

The edge states in 2D TIs have perfectly conducting channels like quantum Hall edge 
states. It is therefore expected to show anomalous thermoelectric transport as well. In 
Ref. 



59] we proposed the 2D TI to be a potential candidate to have efficient thermo- 
electric transport at low temperature. Efficiency of thermoelectric energy conversions is 
characterized by the thermoelectric figure of merit ZT . It is defined as ZT — aS 2 T '/ k, 
where T is the temperature, a is the electrical conductivity, S is the Seebcck coefficient, 
and k is the thermal conductivity. It is not easy to achieve large values of ZT, because 
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a, S and k cannot be independently controlled. We note here that some of the TIs such 



as Bii_ x Sb^ 



6l| . Bi2Se3 5l| . and Bi2Te3 [52j are also good thermoelectric materials. 



To calculate thermoelectric figure of merit, we introduce linear response of the electric 
current j and thermal current w to the electric field or thermal gradient in the following 
way: 

j/q \ _ I L Li 
w J y Li L 2 

where q is the charge — e, and [i is the chemical potential. From the transport coefficients 
Lj , various thermal and electric properties are expressed as 



dp, 

dx 

1 dT 
T dx 



(4) 



1 Li 1 L L 2 - L\ L\ 

a = e Lq, b = , K e = , Zl = s , (5) 

eT L Q T L L L 2 — Lf + k l TL 



where K e and Kl are the thermal conductivity from electrons and phonons, respectively. 
We treat kl as a constant. 

3.2.1. Two-dimensional topological insulators 

In Rcf. 5^] , we showed that the edge states in a narrow ribbon of 2D TI give rise to 
large ZT at low temperature. The edge channels of the 2D TIs are perfectly conducting, 
and it can be described by the Landauer formula with the transmission probability being 
unity. The density of states are schematically shown in Fig.jTJa), where A is the energy 
gap. We focus on energies around the edge of the bulk conduction band and neglect the 
valence band. In the calculation, bulk and edge transports are treated independently, 
which is valid within the inelastic scattering length i. Therefore, we can regard the 
inelastic scattering length £ as an effective system size. The transport coefficients for the 
bulk are calculated by the Boltzmann equation. 

When bulk and edge transport coexist, the total transport coefficients L v are given 
by their sum. The Scebcck coefficients have opposite signs for the bulk and the edge 
transport, and that they cancel each other. As a result the bulk and edge thermoelectric 
transports compete each other, suppressing the ZT in total [g^]. The results are shown 
in Fig. [7] Maximum of ZT occurs when fi is around the band edge because of the 
competition. We note that the edge transport coefficients are proportional to the inelastic 
scattering length £, while the bulk transport does not depend on £. 
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At high temperature such as room temperature, i is very short, and the bulk transport 
is dominant in the thermoelectric transport. At low temperature, the edge transport be- 
comes prominent because £ becomes long, and thermoelectric figure of merit is enhanced. 
We note that for the perfect conduction for edge transport, the ribbon should be much 
wider than the penetration depth A of the edge states. The penetration depth A depends 
on the systems, and in some systems such as Bi ultrathin film, it is as short as the lattice 
constant [63 1. 
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Figure 7: (a) Schematic energy bands for the bulk and edge states used in the calculation, (b) shows 
the results. We employ physical quantities for Bi2Te3 as follows. The electron effective mass is 0.02m e 
where m e is the free electron mass and the number of carrier pockets in the bulk is c = 6. The mobility 
fi* is assumed to be 2000cm 2 V _1 s~ 1 at T =1.8K , and k l is 0.1 Wm^r 1 . The width of the ribbon 
is lOnm. 

In 3D TI systems, topologically protected ID states are predicted on dislocations of 
the crystal, depending on the Zi topological numbers These ID states are perfectly 



conducting, and in Rcf. [65| we proposed that they lead to enhanced thermoelectric figure 
of merit at low temperature, especially for densely distributed dislocations in the crystal. 
The mechanism is similar to that for the edge states of 2D TIs. 



12 



3.2.2. Three- dimensional topological insulators 

Wc discuss here thermoelectric properties of topological surface states in 3D TIs 
However, the surface states in 3D TIs are diffusive, unlike the edge states in 2D TI. 
That reduces the surface transport, and leads to low thermoelectric figure of merit ZT 
due to the competition with the bulk transport. Wc consider a thin slab 3D TI, and 
calculate transport coefficients along the slab for the surface states and the bulk states 
separately by the Boltzmann equation. The phonon transport is neglected. An example 
of the result for ZT is shown in Fig [5] (b). The result shows a competition between the 
surface and bulk states in the thermoelectric transport. The surface transport dominates 
thermoelectricity at low temperatures, and as the temperature is increased, the maximum 
value of the ZT approaches the band edge. 
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Figure 8: (a) A schematic figure of the 3D TI band structure with gapless surface states, and (b) ZT as 
a function of the chemical potential. 



4. Discussion 

4-.1. Finite-size effect in topological insulators 

Most of the intriguing effects discussed above stem from the (massless) Dirac fermion 
on the surface of TIs. However, in general, bulk electrons also contribute to the physics, 
which may mask the predicted effects. Therefore, to reduce the bulk contribution in 
transport properties is one of most important experimental tasks 2J| . To fabricate thin 



films of TIs seems most promising route to overcome this difficulty. However, one should 
keep in mind that when two surface states are sufficiently close to each other, namely 
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Figure 9: Plot of the finite-size induced gap A between the surface states for Bi2Se3 as a function of the 
width L. 



within the decay length of the surface states, these states are overlapped and a gap opens 
in the surface states. [46 . 
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In Ref.[46[, it has been clarified how the surface states in the strong TI Bi2Se3 are 
influenced by finite size effects using the effective Hamiltonian for Bi2Se3. It is found 
that the surface-states in Bi2Se3 display a remarkable robustness towards decreasing the 
width L down to a few nm, thus ensuring that the topological surface states remain 
intact, and that the gapping due to the hybridization of the surface states features an 
oscillating exponential decay as a function of L in Bi2Se3 as shown in Fig. [5) 

To reduce the bulk contribution, gate voltage may be useful, In addition to carrier 



50(, the formation of one-dimensional 



tuning, gate voltage may form ID states. In Rcf. 
channels on the topological surface under the gate electrode has been predicted. This can 
be understood by regarding a momentum in the perpendicular direction as a "mass" . The 
energy dispersion of these channels is almost linear in the momentum, and its velocity 
and sign sensitively depend on the strength of the gate voltage. As a result, the local 
density of states near the gated region has an asymmetric structure with respect to zero 
energy. In the presence of the electron-electron interaction, the correlation effect can be 
tuned by the gate voltage. 
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4-2. Band structure engineering of interface states of topological insulators 

So far we have seen various novel properties of spin transport in topological insulator 
surfaces. The surface band structure typically has a single Dirac cone, with its velocity 
determined by the material itself. The band structure can be designed further, by making 



interfaces between two TIs, as we see in the following [59(. We reconsider the refraction 
problem at the junction between two TIs in Section 2. When the velocities of the two TIs 
have opposite signs. Then both reflection and transmission for the surface transport is 
prohibited for normal incidence, because they violate spin conservation (see Fig. llOf a)). 
Our answer to this paradox is the following. Here we show existence of g apless states at 
the interface between the two TIs (the purple region in Fig. ITUT a)) |59(. The normally 
incident wave goes along the surface of one TI, then flows into the interface. 



(a) 



(b! 




Figure 10: (a) Schematic figure of refraction on the junction between the two TIs (Til (red) and TI2 
(blue)), whose velocities have different signs. The purple region represents the interface between the 
two TIs. (b) For the interface (purple), we consider mirror symmetry with respect to the yz plane, 
and some kind of rotations around the z axis. (c)(d) Dispersion on the interface between the two TIs 
with different signs of velocities. In (c), no rotational symmetry is imposed while in (d) the continuous 
rotational symmetry around the z axis is imposed. 

These interface states arise from hybridization between the two surface states from 

the two TIs. To show the existence of gapless interface states, we first write down the 
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effective Hamiltonian at the interface from the two Dirac cones with hybridization: 

H=\ 1 ), ff 1 =«i(<TXk), ) H 2 = -v 2 (axk) z (6) 

Here #1(2) is the effective surface Hamiltonian for the surface of Til (TI2) at the inter- 
face, and V is the hybridization. For simplicity, we retain only the lowest order in k. We 
henceforth impose the mirror symmetry with respect to the yz plane, M yz , which re- 
stricts the form of the hybridization matrix V . A direct calculation from the Hamiltonian 
(|6]), we can see that for v\v 2 > the interface states are gapped by the hybridization, 
and for v\v 2 < there are gapless states on the interface. Dispersion of the gapless states 
depends on spatial symmetries: 

(i) No rotational symmetry: the interface states have two Dirac cones (Fig. 1101 (d)). 

(ii) Continuous rotational symmetry: the gap of the interface states is closed along a 
loop in k space (Fig. [101 (c)). 

(iii) Discrete crystallographic rotational symmetry: the interface states have a collection 
of Dirac cones, whose number depends on details of the rotational symmetry. For 
C3 symmetry there are six Dirac cones. 

The existence of gapless interface states can also be shown in the following way. Each 
TI with mirror symmetry Ai yz is characterized by the mirror Chern number [53l |. In 
Til and TI2, the two surface Dirac cones have opposite velocities, and it corresponds 
to different mirror Chern numbers; the Til has = — 1, and the TI2 has n T{ = 1. 
This difference by two gives rise to the two Dirac cones on the mirror plane for the 
interface states. These interface states do not come from the Z 2 topological number, but 
come from the mirror Chern number, and are protected by the mirror symmetry M yz . 
Therefore it is a new class of the surface states originated from the mirror symmetry, 
similar to the topological crystalline insulators 
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Both in the present work of 
56|, l57[, the interface/surface 



interface states and in the topological crystalline insulators 
states have even number of Dirac cones, and they are protected by mirror symmetry. In 
the present work Dirac cones appear on an interface between two i^-nontrivial TIs, while 
in the topological crystalline insulator they appear on a surface of a ^-trivial insulator. 
We note that in a related work 581 on an interface between two TIs with the same 



signs of velocities, a surface-state dispersion is reported to have an infinite slope (ve- 
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locity), which is called tachyons. We think that such an infinite slope cannot occur in 
TI interfaces, since the velocity given by v = = |v| "P) (|^): eigenstate) cannot 
diverge even for surface states. The whole eigenvalue problem (i.e. the Schrodingcr 
equation plus decaying boundary conditions at z — > oo) is definitely Hermitian. This 
hermiticity forbids appearance of such tachyon-like modes in the TI junction problems. 

5. Summary 

In this paper, we have investigated spintronics and spincaloritronics in topological 
insulators. We showed spintronics effects in 2D topological insulator junctions and 3D 
topological insulators coupled to ferromagnets. We also investigated spin polarization on 
the surface of a topological insulator induced by a circularly polarized light and refraction 
at the junction between two topological insulators. As for spincaloritronics effects, we 
showed transverse magnetic heat transport and thermoelectric transport in topological 
insulators. Finite-size effect in topological insulators and band structure engineering of 
interface states of topological insulators are also discussed. 
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